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Abstract 



We discuss the electromagnetic-penguin-dominated radiative B decays 
B ^ Xs + 'y, K*^{o) + and 5^ ^ + 7 in the context of 

the standard model (SM) and their Cabibbo-Kobayashi-Maskawa (CKM)- 
suppressed counterparts, B — > Xd + 7, -B^ — + 7, — {p^,uj) + 7, 
and Bs K*^ + 7, using QCD sum rules for the exclusive decays. The 
importance of these decays in determining the parameters of the CKM matrix 
[|^] is emphasized. The semileptonic decays B — > Xsi~^i~ are also discussed in 
the context of the SM and their role in determining the Wilson coefficients of 
the effective theory is stressed. Comparison with the existing measurements 
are made and SM-based predictions for a large number of rare B decays are 
presented. 
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1 Estimates of B{B ^ + 7) and \Vts\ in the 
Standard Model 

The Standard Model (SM) of particle physics does not admit Flavour-changing 
neutral-current (FCNC) transitions in the Born approximation. However, 
they are induced through the exchange of bosons in loop diagrams. The 
short- distance contribution in rare decays is dominated by the (virtual) top 
quark contribution. Hence the decay characteristics provide quantitative 
information on the top quark mass and the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements Vu, i = d, s,b We shall discuss representative 
examples from several such transitions involving B decays, starting with the 
decay i? — > + 7, which has been measured by CLEO [§. This was pre- 
ceded by the measurement of the exclusive decay B — > K* + 7 by the same 
collaboration ^j. The present measurements give [Q: 

i3(5^X, + 7) = (2.32 ± 0.57 ± 0.35) X 10~\ (1) 
BiB^K^ + ^f) = (4.2±0.8±0.6) X 10~\ (2) 

yielding an exclusive-to-inclusive ratio: 

«- = ^^|^5J^^ = (l"±6-«)%. (3) 

These decay rates determine the ratio of the CKM matrix elements |VJs|/|V^6| 
and the quantity Rk* provides information on the decay form factor in B ^ 
K* + 7. In what follows we take up these points briefly. 

The leading contribution tob —* s+7 arises at one-loop from the so-called 
penguin diagrams. With the help of the unitarity of the CKM matrix, the 
decay matrix element in the lowest order can be written as: 

Mib-^s = ^^\t{F2{xt)-F2{x,))q^t^-sa^,{m,R + m,L)b . (4) 
V2 

where Xi = mf/niyy, and and are, respectively, the photon four- 
momentum and polarization vector. The GIM mechanism |Q is manifest 
in this amplitude and the CKM-matrix element dependence is factorized in 
Xf = VtbV^*. The (modified) Inami-Lim function F2{xi) derived from the 
(1-loop) penguin diagrams is given by ||]: 

F^{x) = [6x(3a;-2)logx-(x-l)(8x2 + 5x-7)l . (5) 
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The measurement of the branching ratio for 5 — >■ + 7 can be readily 
interpreted in terms of the CKM-matrix element product At/|Vcfe| or equiv- 
alently | Vj^l/I For a quantitative determination of |V^s|/|Kfe|, however, 
QCD radiative corrections have to be computed and the contribution of the 
so-called long-distance effects estimated. 

The appropriate framework to incorporate QCD corrections is that of 
an effective theory obtained by integrating out the heavy degrees of freedom, 
which in the present context are the top quark and bosons. The operator 
basis depends on the underlying theory and for the SM one has (keeping 
operators up to dimension 6), 

Heffib ^ . + 7) = -^v,:v,,J2c^i^^)o^i^^), (6) 

where the operator basis, the "matching conditions" Ci{mw), and the solu- 
tions of the renormalization group equations Ci{fi) can be seen in ref. [0]. 
The perturbative QCD corrections to the decay rate r{B + 7) have 

two distinct contributions: 

• Corrections to the Wilson coefficients Ci{fi), calculated with the help of 
the renormalization group equation, whose solution requires the knowl- 
edge of the anomalous dimension matrix in a given order in a^. 

• Corrections to the matrix elements of the operators Oi entering through 
the effective Hamiltonian at the scale n = 0{mh). 

The anomalous dimension matrix is needed in order to sum up large log- 
arithms, i.e., terms like a" (mw) log™ (mfe/M), where M = rrit or mw and 
m < n (with n = 0, 1, 2, ...). At present only the leading logarithmic correc- 
tions (m = n) have been calculated systematically and checked by several 
independent groups in the complete basis given in Eq. (|^) First calcula- 
tions of the NLO corrections to the anomalous dimension matrix have been 
recently reported by Misiak and are found to be small. Next-to-leading 
order corrections to the matrix elements are now available completely. They 
are of two kinds: 

• QCD Bremsstrahlung corrections b ^ s'-f + g, which are needed both 
to cancel the infrared divergences in the decay rate for B ^ Xg + j 
and in obtaining a non-trivial QCD contribution to the photon energy 
spectrum in the inclusive decay B Xs + 
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Next-to-leading order virtual corrections to the matrix elements in the 
decay 6 — > s -|- 7. 



The Bremsstrahlung corrections were calculated in [[I^, |ri|] in the truncated 
basis and last year also in the complete operator basis |T2| , |13[. The higher 
order matching conditions, i.e., Ci{mw), are known up to the desired accu- 
racy, i.e., up to 0{as{Miy)) terms WM. The next-to- leading order virtual 



corrections have also been calculated |jT5[. They reduce the scale-dependence 



of the inclusive decay width. The branching ratio B{B + 7) can be 

expressed in terms of the semileptonic decay branching ratio 

B{B ^ X,7) = [£CB_^2±^]t/. siB Xeu,) , (7) 

J- SL 

where the leading-order QCD corrected expression for Tsl can be seen in [0]. 
The leading order (1/mb) power corrections in the heavy quark expansion are 
identical in the inclusive decay rates for B Xg + 'j and B Xiui, entering 
in the numerator and denominator in the square bracket, respectively, and 
hence drop out. 

In Ref. [|^, the present theoretical errors on the branching ratio B{B — > 
Xg'y) are discussed, yielding: 

B{B ^ X, + 7) = (3.20 ± 0.30 ± 0.38 ± 0.32) x 10"^ (8) 

where the first error comes from the combined effect of Arrit and Afi (the 
scale dependence), the second error arises from the extrinsic sources (such as 
A(mfe), A{BRsl)), and the third error is an estimate (±10%) of the NLO 
anomalous dimension piece in C^-^ , the coefficient of the magnetic moment 
operator. Combining the theoretical errors in quadrature gives [0]: 

i3(5 ^ X, + 7) = (3.20 ± 0.58) X 10-^ (9) 

which is compatible with the present measurement B{B — Xs + 7) = (2.32 ± 
0.67) X 10~^ p|. Expressed in terms of the CKM matrix element ratio, one 

gets 

= 0.85±0.12(expt) ±0.10(th), (10) 



IK 



cb\ 



which is within errors consistent with unity, as expected from the unitarity 
of the CKM matrix. 
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2 Inclusive radiative decays B + 7 

The theoretical interest in studying the (CKM- suppressed) inclusive radia- 
tive decays B — > Xd + 7 lies in the first place in the possibility of deter- 
mining the parameters of the CKM matrix. We shall use the Wolfenstein 
parametrization [|l^], in which case the matrix is determined in terms of the 



four parameters A, X = sin6'c, p and r]. The quantity of interest in the de- 
cays B + 7 is the end-point photon energy spectrum, which has to be 
measured requiring that the hadronic system recoiling against the photon 
does not contain strange hadrons to suppress the large- i?^ photons from the 
decay i? — > + 7. Assuming that this is feasible, one can determine from 
the ratio of the decay rates B{B ^ X^ + l)/B{B ^ X, + 7) the CKM- 
Wolfenstein parameters p and rj. This measurement was first proposed in 



[|TT|], where the photon energy spectra were also worked out. 

In close analogy with the B — >■ Xs-l-7 case discussed earlier, the complete 
set of dimension-6 operators relevant for the processes b ^ d'y and b d''yg 
can be written as: 

- ^) = 6 E (11) 

where = Vjb V*^ for j = t, c, u. The operators Oj, j = 1, 2, have implicit in 
them CKM factors. In the Wolfenstein parametrization |T^, one can express 
these factors as : 

^u = AX'{p-t7]), ^c = ~AX\ et = -e«-ec. (12) 

We note that all three CKM-angle-dependent quantities C,j are of the same 
order of magnitude, O(A^). This aspect can be taken into account by defining 
the operators Oi and O2 entering in TLeff{b d) as follows |11|]: 



01 = --^{cLfil^bLa){,dLal,iCLp) - -^{uLfil^bLa){,dLa'l i^Uip) , 

■St Kt 

02 = ~{cLal^hLa){dLpli,CLp) - ^{uLal'^bia) {dLfSlf^ULfs) , (13) 

with the rest of the operators [Oj] j = 3. ..8) defined like their counterparts 
Oj in Tieffip — > s), with the obvious replacement s —>■ d. With this choice, the 
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matching conditions Cj{mw) and the solutions of the RG equations yielding 
Cj(/i) become identical for the two operator bases Oj and Oj. The essential 
difference between r{B + 7) and r{B —>■ + 7) lies in the matrix 

elements of the first two operators Oi and O2 (in Tieffib s)) and Oi and 
O2 (in Tieffib —>■ d)). The branching ratio B{B —>■ + 7) in the SM can be 
written as: 

BiB ^ X, + 7) = D,X'{{1 - pf + r^^ - (1 - p)D2 - vD-s + D,}, (14) 

where the functions Di depend on the parameters mt,mb,mc, p, as well as 
the others we discussed in the context of B{B — + 7). These functions 



were first calculated in |]TT| in the leading logarithmic approximation. Re- 



cently, these estimates have been improved in |T^, making use of the NLO 



calculations in [|T5|. To get an estimate of the inclusive branching ratio, the 
CKM parameters p and r] have to be constrained from the unitarity fits. 
Present data and theory restrict them to lie in the following range (at 95% 
C.L.) 0: 

0.20 < ri< 0.52, 
-0.35 < p< 0.35 , (15) 

which, on using the current lower bound from LEP on the B^ - B^ mass 
difference AM, > 9.2 (ps)(-^) [gUf, is reduced to -0.25 < P < 0.35 using 
= 1-1, where C,s is the 5'f/(3)-breaking parameter = fs^BBs/ fB^BBa- 
The preferred CKM-fit values are {p,rj) = (0.05,0.36), for which one gets 



S(S ^ Xrf + 7) = 1.63 X 10-^ (16) 

whereas B{B ^ + 7) = 8.0 x 10"^ and 2.8 x 10~^ for the choice p = 
0.35, T] = 0.40 and p = —77 = —0.25, respectively. In conclusion, we note that 
the functional dependence of B{B X^ + 7) on the Wolfenstein parameters 
(p, 77) is mathematically different than that of AM,. However, qualitatively 
they are very similar. From the experimental point of view, the situation 
p < is favourable for both the measurements as in this case one expects 
(relatively) smaller values for AM, and larger values for the branching ratio 
B{B Xd + 7), as compared to the p > case which would yield larger 
AM, and smaller B{B ^ X^ + 7). 
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2.1 B{B V + j) and constraints on the CKM param- 
eters 

Exclusive radiative B decays B —>■ V + 'j, with V = K*,p,uj, are also po- 
tentially very interesting from the point of view of determining the CKM 
parameters The extraction of these parameters would, however, in- 

volve a trustworthy estimate of the SD- and LD-contributions in the decay 
amplitudes. 

The SD-contribution in the exclusive decays {B^, B^) {K*^, K*^) +7, 
{B^,B^) — > {p^,p^) + 7, — i> + 7 and the corresponding Bg decays, 
-Bs ^ + 7, and Bs K*^ + 7, involve the magnetic moment operator 
Oj and the related one obtained by the obvious change s — > c?, O7. The 
transition form factors governing the radiative B decays i? ^ F + 7 can be 
generically defined as: 

(y, \\^-i,a,,q''h\B) = 2e^.p.eWp^p^Fj^^^(0). (17) 

Here is a vector meson with the polarization vector e^'^\ V = p,uj,K* or 
0; -B is a generic 5-meson B^, B^ or Bg, and ip stands for the field of a light 
u,d or s quark. The vectors Pb, Pv and q = Pb — Pv correspond to the 4- 
momenta of the initial 5-meson and the outgoing vector meson and photon, 
respectively. In (^) the QCD renormalization of the ipcr^u^fb operator is 
implied. Keeping only the SD-contribution leads to obvious relations among 
the exclusive decay rates, exemplified here by the decay rates for Bd) — *■ 
p + 7 and Bd) ^K*+r. 
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■ u,d — ^uA 



ts\ 



r((5±,i?0)-^(K*±,K*o)+7) 

(18) 

where ^u,d is a phase-space factor which in all cases is close to 1 and Hi = 
[Fs{Bi p'y)/Fs{Bi K*'y)Y. The transition form factors Fs are model 
dependent. Estimates of Fs in the QCD sum rule approach in the normal- 
ization of Eq. (^) range between Fs{B K*'y) = 0.31 (Narison in p2|) to 
Fs{B K*-i) = 0.37 (Ball in ||), with a typical error of ±15%, and hence 
are all consistent with each other. This, for example, gives Rk* = 0.16±0.05, 



using the result from which is in good agreement with data. The ra- 



tios of the form factors, i.e. Kj, should therefore be reliably calculable as 
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they depend essentially only on the SU(3)-breaking effects which have been 
estimated 
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The LD-amplitudes in radiative B decays from the light quark intermedi- 
ate states necessarily involve other CKM matrix elements. Hence, the simple 
factorization of the decay rates in terms of the CKM factors involving \Vtd\ 
and \Vts\ no longer holds thereby invalidating the relation ([18|) given above. 
In the decays B ^ V + j they are induced by the matrix elements of the 
four-Fermion operators Oi and O2 (likewise Oi and O2). Estimates of these 
contributions require non-perturbative methods. This problem has been in- 
vestigated in ^ using a technique |^ which treats the photon emission 
from the light quarks in a theoretically consistent and model-independent 
way. This has been combined with the light-cone QCD sum rule approach 
to calculate both the SD and LD — parity conserving and parity violating 
— amplitudes in the decays {B^,B^) {p^,pluj) +7. To illustrate this, we 
concentrate on the B^ decays, B^ — > and take up the neutral B decays 

B^ — > p{uj) -|- 7 at the end. The LD-amplitude of the four-Fermion operators 
Oi, O2 is dominated by the contribution of the weak annihilation of valence 
quarks in the B meson and it is color-allowed for the decays of charged 5^ 
mesons. Using factorization, the LD-amplitude in the decay Bu 
can be written in terms of the form factors and F2 , 



long 



uuaf3 



2F^(g2)}(.19) 



Again, one has to invoke a model to calculate the form factors. Estimates 
from the light-cone QCD sum rules give ||24|] : 



F^/Fs = 0.0125 ±0.0010; 



0.0155 ±0.0010 



(20) 



where the errors correspond to the variation of the Borel parameter in the 
QCD sum rules. Including other possible uncertainties, one expects an ac- 
curacy of the ratios in ( PO] ) of order 20%. The parity-conserving and parity- 
violating amplitudes turn out to be numerically close to each other in the 
QCD sum rule approach, F(" ~ = Fl, hence the ratio of the LD- and the 
SD- contributions reduces to a number \UM 



•^long / 



short 



R 



L/S 



Vt,v:, 



(21) 
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Using C2 = 1.10, Ci = -0.235, Cf = -0.306 from Ref. [0 (corresponding 
to the scale /z = 5 GeV) gives: 



L 



i?" ^ = ' '•^ ■ ± = -0.30 ± 0.07 , (22) 

which is not small. To get a ball-park estimate of the ratio Aiongl Ashort^ we 
take the central value from the CKM fits, yielding |Vufe|/|VJd| — 0.33 |]lj 



Along/ AsHortr^" ^ = \R^L/S^' 'l^f^ ^ 10% . (23) 



^tdVtb 

Thus, the CKM factors suppress the LD-contributions. 

The analogous LD-contributions to the neutral B decays p7 and 

wf are expected to be much smaller. The corresponding form factors 
for the decays B^ — >■ p°(a;)7 are obtained from the ones for the decay B^ — >■ 
p^7 discussed above by the replacement of the light quark charg 
which gives the factor —1/2; in addition, and more importantly, the LD- 
contribution to the neutral B decays is colour-suppressed, which reflects itself 
through the replacement of the factor ai by 02. This yields for the ratio 



R 



L/S _ 6^02 



eua^ 



-0.13 ±0.05, (24) 



where the numbers are based on using 02/01 = 0.27 ± 0.10 This would 
then yield R^j^'^ — R^/'s^"' = 0.05, which in turn gives 

< 0.02. (25) 

•^short 

This, as well as the estimate in eq. P3, should be taken only as indicative 



in view of the approximations made in [p3| , p4| . That the LD-effects remain 
small in B^ p7 has been supported in a recent analysis based on the soft- 
scattering of on-shell hadronic decay products B^ p^p^ p'-f [^, though 
this paper estimates them somewhat higher (between 4 — 8%). 

Restricting to the colour-allowed LD-contributions, the relations, which 
obtains ignoring such contributions (and isospin invariance), 

r(fi=^ ^ p±7) = 2 r(5° ^ p%) = 2 r(5° ^ u-f) , (26) 
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get modified to 



1 I pB±^p±7 KfeKd 



21(5° p7) 21(5° ^ 

= 1 + 2 . _ ^ + + • (27) 

where i?L/5 = Rf/^"^^ ■ The ratio r(5± ^ p±7)/2r(50 ^ p7)(= r(S± ^ 
p^7)/2r(i?° — i> 077)) is shown in Fig. |l] as a function of the parameter p, 
with 7] = 0.2, 0.3 and 0.4. This suggests that a measurement of this ratio 
would constrain the Wolfenstein parameters (p, t]) , with the dependence on 
p more marked than on rj. In particular, a negative value of p leads to a 
constructive interference in Bu —>■ P7 decays, while large positive values of p 
give a destructive interference. 

The ratio of the CKM-suppressed and CKM-allowed decay rates for charged 
B mesons gets modified due to the LD contributions. Following earlier dis- 
cussion, we ignore the LD-contributions in r(i? — > K*'-/). The ratio of the 
decay rates in question can therefore be written as: 

r(i?±-.p±7) ,2rn , .1 



r(5± ^ K*^-f) 



X 



Using the central value from the estimates of the ratio of the form factors 
squared k„ = 0.59 ± 0.08 [^, we show the ratio (pSD in Fig. || as a function 



of p for 7] = 0.2, 0.3, and 0.4. It is seen that the dependence of this ratio is 
rather weak on rj but it depends on p rather sensitively. The effect of the 
LD-contributions is modest but not negligible, introducing an uncertainty 
comparable to the ~ 15% uncertainty in the overall normalization due to the 
S'[/(3)-breaking effects in the quantity Ku- 

Neutral 5-meson radiative decays are less-prone to the LD-effects, as 
argued above, and hence one expects that to a good approximation (say, 
better than 10%) the ratio of the decay rates for neutral B meson obtained 
in the approximation of SD-dominance remains valid [^: 



r{B ^ K*-f) 



n,X'[{l-py + r^'] , (29) 
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Figure 1: Ratio of the neutral and charged -B-decay rates r(_B='= 
p^7)/2r(i?° — i> fTf) as a function of the Wolfenstein parameter p, with 
T] = 0.2 (short-dashed curve), t] = 0.3 (sohd curve), and r] = 0.4 (long-dashed 
curve). (Figure taken from P^.) 



where this relation holds for each of the two decay modes separately. 

Finally, combining the estimates for the LD- and SD-form factors in |^ 
and ||2T| , respectively, and restricting the Wolfenstein parameters in the range 
-0.25 < p< 0.35 and 0.2 < r/ < 0.4, as suggested by the CKM-fits 0, we 
give the following ranges for the absolute branching ratios: 

B{B^ ^ p^-f) = (1.5± 1.1) X 10-^ , 
B{B^ p-i) - B{B^ CU7) = (0.65 ± 0.35) x 10"^ , (30) 

where we have used the experimental value for the branching ratio B{B — >■ 
^* + 7) 0) adding the errors in quadrature. The large error reflects the 
poor knowledge of the CKM matrix elements and hence experimental deter- 
mination of these branching ratios will put rather stringent constraints on 
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Figure 2: Ratio of the CKM-suppressed and CKM-allowed radiative B- 
decay rates T{Bu — >■ pj)/T{B — > K*j) (with B = B^ or B^) as a function of 
the Wolfenstein parameter p, a) with r] = 0.2 (short-dashed curve), t] = 0.3 
(sohd curve), and r] = 0.4 (long-dashed curve). (Figure taken from P^ .) 



the Wolfenstein parameter p. 

In addition to studying the radiative penguin decays of the B^ and 5^ 
mesons discussed above, hadron machines such as HERA-B will be in a po- 
sition to study the corresponding decays of the B^ meson and Af, baryon, 
such as — i> + 7 and Af, — A + 7, which have not been measured so far. 
We list below the branching ratios in a number of interesting decay modes 



calculated in the QCD sum rule approach in |21 



B{Bs 07) = B{Bd K*-f) = (4.2 ± 2.0) x 10 
- (0,36 ±0,14)1 1''"!- 



-5 



=^ i3(fi, ^ ir*7) = (0.75 ± 0.5) X lO"*' . (31) 
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The estimated branching ratios in a number of inclusive and exclusive 
radiative B decay modes are given in Table |I|, where we have also listed the 
branching ratios for Bg 77 and 5^ — 77. 

2.2 Inclusive rare decays B Xgt^l^ in the SM 

The decays B —>■ Xgi'^i' , with i = e, fi, r, provide a more sensitive search 
strategy for finding new physics in rare B decays than for example the decay 
B Xsj , which constrains the magnitude of Cj-^ . The sign of Cj-^ , which 
depends on the underlying physics, is not determined by the measurement 
of B{B — i> + 7). This sign, which in our convention is negative in the 
SM, is in general model dependent. It is known (see for example [^) that 
in supersymmetric (SUSY) models, both the negative and positive signs are 
allowed as one scans over the allowed SUSY parameter space. We recall 
that for low dilepton masses, the differential decay rate for B — > Xsi~^i~ is 
dominated by the contribution of the virtual photon to the charged lepton 
pair, which in turn depends on the effective Wilson coefficient C^'^ . However, 
as is well known, the B Xsi~^i~ amplitude in the standard model has two 
additional terms, arising from the two FCNC four-Fermi operators, which 
are not constrained by the 5 ^ + 7 data. Calling their coefficients Cg 
and Cio, it has been argued in that the signs and magnitudes of all three 
coefficients C^-^, Cg and Cio can, in principle, be determined from the decays 
5 X, + 7 and 5 ^ XJ+i- . 

The SM-based rates for the decay b — > si^i~ , calculated in the free quark 
decay approximation, have been known in the LO approximation for some 



time [^. The LO calculations have the unpleasant feature that the decay 
distributions and rates are scheme-dependent. The required NLO calcula- 
tion is in the meanwhile available, which reduces the scheme-dependence of 
the LO effects in these decays In addition, long-distance (LD) effects, 



which are expected to be very important in the decay B Xsi^i~ , have 
also been estimated from data on the assumption that they arise dominantly 
due to the charmonium resonances J/ip and ip' through the decay chains 
B XsJ/ip{ip', ...) Xsi^i~ . Likewise, the leading {l/rrib^) power correc- 
tions to the partonic decay rate and the dilepton invariant mass distribution 

^This also holds for a large class of models such as MSSM and the two-Higgs doublet 
models but not for all SM-extensions. In LR symmetric models, for example, there are 
additional FCNC four- Fermi operators involved 
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have been calculated with the help of the operator product expansion in the 
effective heavy quark theory [^]. The results of |31| have, however, not 
been confirmed in a recent independent calculation , which finds that the 
power corrections in the branching ratio B{B Xsi^i~) are small (typically 
— 1.5%). The corrections in the dilepton mass spectrum and the FB asymme- 
try are also small over a good part of this spectrum. However, the end-point 
dilepton invariant mass spectrum is not calculable in the heavy quark expan- 
sion and will have to be modeled. Non-perturbative effects in B ^ Xsi~^i~ 
have been estimated using the Fermi motion model in These effects 



are found to be small except for the end-point dilepton mass spectrum where 
they change the underlying parton model distributions significantly and have 
to be taken into account in the analysis of data [B^ . 

The amplitude for B Xs£^i~ is calculated in the effective theory ap- 
proach, which we have discussed earher, by extending the operator basis of 
the effective Hamiltonian introduced in Eq. (0): 



Heffib^ s + r,b^ s + i+r) 

n,ff{b^s + ^)-^v*v,,[c,o. 



10>-^10j ! 



(32) 



where the two additional operators are: 



0,0 



a 
a 



SaYPlbJ 



(33) 



The analytic expressions for d^irriw) and CiQ{mw) can be seen in 



and will not be given here. We recall that the coefficient Cg in LO is scheme- 
dependent. However, this is compensated by an additional scheme- dependent 
part in the (one loop) matrix element of Og. We call the sum C^^ , which is 
scheme-independent and enters in the physical decay amplitude given below, 



M{b 



s + 



CfsYPLb£i,£ + CiosYPLbi^^.7,1 



2Cf-sia^,%{m,PR + m,P,.)6£7^£ 



(34) 
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with 

Cf{s)^C,r]is)+Yis). (35) 

The function Y{s) is the one-loop matrix element of Og and can be seen in 
literature 0]. A useful quantity is the differential FB asymmetry in the 
c.m.s. of the dilepton defined in refs. [p5[] : 



dAis) _ /■! dB /--i dB 

(36) 



ds Jo dz Jo dz 

where z = cos 6, with 6 being the angle between the lepton and the 
6-quark. This can be expressed as: 

^ = -i3.,|^^t.^(S)Cio[5K(Cf (5)) + 2C,^^^(l + m^)].(37) 

The Wilson coefficients C^'^^ , C^^^ and Cio appearing in the above equation 
and the dilepton spectrum (see, for example |^) can be determined from 
data by solving the partial branching ratio i3(As) and partial FB asymmetry 
A{/S.s), where As defines an interval in the dilepton invariant mass PH| . 
There are other quantities which one can measure in the decays B — >■ 
to disentangle the underlying dynamics. We mention here the lon- 
gitudinal polarization of the lepton in i? ^ Xs£^£~ , in particular in i? — *^ 
XsT^T^ , proposed by Hewett In a recent paper, Kriiger and Sehgal 



have stressed that complementary information is contained in the two or- 
thogonal components of polarization {Pt-, the component in the decay plane, 
and P/v, the component normal to the decay plane), both of which are pro- 
portional to m^/mfc, and therefore significant for the r+r" channel. A third 
quantity, called energy asymmetry, proposed by Cho, Misiak and Wyler [P^], 
defined as 

. ^ N{E,->E,^)~N{E,+ >E,-) 

N{E,- > E,+ ) + N{E,+ > E,-) ' ^ 

where N{Ei- > E£+) denotes the number of lepton pairs where is more 
energetic than i~ in the i?-rest frame, is, however, not an independent mea- 
sure, as it is directly proportional to the FB asymmetry discussed above. 
The relation is ||3^ : 

''Ais)=BxA. (39) 
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This is easy to notice if one writes the Mandelstam variable u{s) in the 
dilepton cm. and the S-hadron rest systems. 

Next, we discuss the effects of LD contributions in the processes B —>■ 
Xsi~^i~. Note that the LD contributions due to the vector mesons such as 
J/ip and ip', as well as the continuum cc contribution already discussed, ap- 
pear as an effective {sL'~^^hL){i'~^^€) interaction term only, i.e. in the operator 
O9. This implies that the LD-contributions should change Cg effectively, C7 
as discussed earlier is dominated by the SD-contribution, and Cio has no 
LD-contribution. In accordance with this, the function Y{s) is replaced by, 

Y{s) ^ Y'{s) = Y{s) + Fres(5), (40) 



where Iresl-s) is given as ||35| 



Yresi-s) = -,k (3Ci + C, + SO, + C, + 3C, + C,) Yl ,5^^^^/^^^^^ 

(41) 

where k is a fudge factor, which appears due to the inadequacy of the fac- 
torization framework in describing data on i? ^ J/%l)Xs. The long-distance 
effects lead to significant interference effects in the dilepton invariant mass 
distribution and the FB asymmetry in i? — > Xsl'^l~ shown in Figs. ^ and 
respectively. This can be used to test the SM, as the signs of the Wilson 
coefficients in general are model dependent. For further discussions we refer 
to Ref. where also theoretical dispersion on the decay distributions due 
to various input parameters is worked out. Taking into account the spread 
in the values of the input parameters, /i. A, m^, and Bsl discussed in the 
previous section in the context of B{B — + 7), we estimate the following 
branching ratios for the SD-piece only (i.e., from the intermediate top quark 
contribution only) |[3^ : 

B{B^X,e+e-) = (8.4 ± 2.3) x 10^^ 
B{B ^ XsjJ^fi-) = (5.7± 1.2) X 10^^ 
B{B^XsT+T-) = (2.6±0.5) X 10"^ (42) 

where theoretical errors and the error on Bsl have been added in quadrature. 
The present experimental limit for the inclusive branching ratio in B ^ 
Xsi~^£~ is actually still the one set by the UAl collaboration some time ago 
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Figure 3: Dilepton invariant mass distribution in B —>■ Xsi~^i~ in the SM 
with the next-to-leading order QCD corrections and non-perturbative effects 
calculated in the Fermi motion model (solid curve), and including the LD- 
contributions (dashed curve). The model parameters {pF,^q) are indicated 
in the figure. Note that the height of the J/ip peak is suppressed due to the 
linear scale. (Figure taken from [B^O 



39|] , namely 13{B Xgfi^fi ) > 5.0 x 10 . As far as we know, there are no 



interesting limits on the other two modes, involving Xse~^e~ and X^r+r^. 

It is obvious from Fig. ^ that only in the dilepton mass region far away 
from the resonances is there a hope of extracting the Wilson coefficients 
governing the short-distance physics. The region below the J/ip resonance is 
well suited for that purpose as the dilepton invariant mass distribution there 
is dominated by the SD-piece. Including the LD-contributions, following 
branching ratio has been estimated for the dilepton mass range 0.2 < s < 0.36 
in H: 

}3{B = (1-3 ± 0.3) X 10"^ (43) 
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Figure 4: Normalized FB asymmetry m.B Xst^i" in the SM as a function 
of the dilepton invariant mass calculated using the next-to-leading order QCD 
correction and the Fermi motion effects (solid curve), and including the LD- 
contributions (dashed curve). The Fermi motion parameters are indicated in 
the figure. (Figure taken from p2[.) 



with B{B XsC^e^) ~ B{B Xs/i+yU"). The FB-asymmetry is estimated 
to be in the range 10% - 27%, as can be seen in Fig. ^ 

The experimental limits on the decay rates of the exclusive decays B — >■ 
{K, K*)i+t 



^ 0, while arguably closer to the SM-based estimates, can 
only be interpreted in specific models of form factors, which hinders some- 
what their transcription in terms of the information on the underlying Wilson 
coefficients. Using the exclusive-to-inclusive ratios 



and 



Rku = T{B Ki+r)/T{B ^ Xs 



Rkhi = T{B K*tr)/T{B ^ X, 



0.07 ±0.02, 



0.27 ±0.0.07, 
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which were estimated in , the results are presented in Table |1|. 

In conclusion, the semileptonic FCNC decays B — > Xsi~^i~ (and also the 
exclusive decays) will provide very precise tests of the SM, as they will deter- 
mine the signs and magnitudes of the three Wilson coefficients, C7, Cg-^ and 
CiQ. This, perhaps, may also reveal physics beyond-the-SM if it is associated 
with not too high a scale. The MSSM model is a good case study where 



measurable deviations from the SM are anticipated and worked out [28 



2.3 Summary and overview of rare B decays in the SM 

The rare B decay mode B — > XgUu, and some of the exclusive channels 
associated with it, have comparatively larger branching ratios. The estimated 
inclusive branching ratio in the SM is [H2[ - P4|| : 



B{B X,vv) = (4.0 ± 1.0) X 10"^ , (44) 

where the main uncertainty in the rates is due to the top quark mass. The 
scale-dependence, which enters indirectly through the top quark mass, has 



been brought under control through the NLL corrections, calculated in 
The corresponding CKM-suppressed decay B — > Xai'u is related by the ratio 
of the CKM matrix element squared [E2| : 



B{B ^ Xduu) 
B{B -> Xsvv) 



td\ 



Wt 



(45) 



ts\ 



Similar relations hold for the ratios of the exclusive decay rates which depend 
additionally on the ratios of the form factors squared, which deviate from 
unity through SU (3)-breaking terms, in close analogy with the exclusive 
radiative decays discussed earlier. These decays are particularly attractive 
probes of the short-distance physics, as the long-distance contributions are 
practically absent in such decays. Hence, relations such as the one in (^) 
provide, in principle, one of the best methods for the determination of the 
CKM matrix element ratio |Vt(i|/|Vis| [^]. From the practical point of view, 
however, these decay modes are rather difficult to measure, in particular at 
the hadron colliders and probably also at the B factories. The best chances 
are in the Z'^-decays at LEP, from where the present best upper limit stems 

B{B Xvv) < 7.7 X 10"^ (46) 
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The estimated branching ratios in a number of inclusive and exclusive decay 
modes are given in Table updating the estimates in 0. 

Further down the entries in Table |I| are listed some two-body rare de- 
cays, such as (-8°, 5°) 77, studied in where only the lowest or- 



der contributions are calculated, i.e., without any QCD corrections, and the 
LD-effects, which could contribute significantly, are neglected. The decays 
B^) — > have been studied in the next-to- leading order QCD in [^] . 
Some of them, in particular, the decays 5° /i and perhaps also the 
radiative decay 77, have a fighting chance to be measured at LHC 

The estimated decay rates, which depend on the pseudoscalar coupling con- 
stant /b^ (for _Bs-decays) and fs^ (for S^-decays), together with the present 
experimental bounds are listed in Table |I|. Since no QCD corrections have 
been included in the rate estimates of {Bg, B^j — > 77, the branching ratios 
are rather uncertain. The constraints on beyond-the-SM physics that will 
eventually follow from these decays are qualitatively similar to the ones that 
(would) follow from the decays 5 — > + 7 and B Xst^l~ ^ which we 
have discussed at length earlier. 



3 Acknowledgements 

I would like to thank Hrachia Asatrian, Vladimir Braun, Christoph Greub 
and Tak Morozumi for helpful discussions. The warm hospitality extended 
by Fernando Ferroni and his collaborators in Rome during the Beauty 96 
workshop is thankfully acknowledged. Finally, I acknowledge the editorial 
assistance of Peter Schlein in the preparation of this manuscript. 



19 



References 

[1] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; 

M. Kobayashi and K. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 

[2] M.S. Alam et al. (CLEO Collaboration), Phys. Rev. Lett. 74 (1995) 
2885. 

[3] R. Ammar et al. (CLEO Collaboration), Phys. Rev. Lett. 71 (1993) 674. 

[4] R. Ammar et al. (CLEO Collaboration), contributed paper to the In- 
ternational Conference on High Energy Physics, Warsaw, 25-31 July 
1996, CLEO CONE 96-05. 

[5] S.L. Glashow, J. Iliopoulos and L. Maiani, Phys. Rev. D2 (1970) 1285. 

[6] T. Inami and C.S. Lim, Prog. Theor. Phys. 65 (1981) 297. 

[7] A. Ah, preprint DESY 96-106 ||hep-ph/ 9606321 ; to appear in the Pro- 
ceedings of the XX International Nathiagali Conference on Physics and 
Contemporary Needs, Bhurban, Pakistan, June 24- July 13, 1995 (Nova 
Science Publishers, New York, 1996). 

[8] M. Ciuchini et al., Phys. Lett. B316 (1993) 127; Nucl. Phys. B415 
(1994) 403; G. Cella et ah, Phys. Lett. B325 (1994) 227; 
M. Misiak, Nucl. Phys. B393 (1993) 23; [E. B439 (1995) 461]. 

[9] M. Misiak, contribution to the International Conference on High Energy 
Physics, Warsaw, 25-31 July 1996. 

[10] A. Ah and C. Greub, Z. Phys. C49 (1991) 431; Phys. Lett. B259 (1991) 
182; Z. Phys. C60 (1993) 433. 

[11] A. Ah and C. Greub , Phys. Lett. B287 (1992) 191. 

[12] A. Ah and C. Greub, Phys. Lett. B361 (1995) 146. 

[13] N. Pott, Phys. Rev. D54 (1996) 938. 



20 



[14 
[15 

[16 

[17; 

[18 
[19 
[20 

[21 
[22 

[23 

[24 
[25 

[26 
[27 



K. Adel and Y.-p. Yao, Phys. Rev. D49 (1994) 4945. 

C. Greub, T. Hurth and D. Wyler, Phys. Lett. B380 (1996) 385; Phys. 
Rev. D54 (1996) 3350. 

L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 

A. Ah and D. London, preprint DESY 96-140, UdeM-GPP-TH-96-45, 
|hep-ph/9607392|| , to appear in the Proc. of QCD Euroconference 96, 



MontpeUier, July 4-12, 1996. 

R.M. Barnett et ah (Particle Data Group), Phys. Rev. D54 (1996) 1. 

A. Ali, H.M. Asatrian, and C. Greub, to be published. 

L. Gibbons, Plenary talk at the International Conference on High Energy 
Physics, Warsaw, ICHEP96 (1996). 

A. Ah, V.M. Braun and H. Simma, Z. Phys. C63 (1994) 437. 

P. BaU, TU-Miinchen Report TUM-T31-43/93 (1993); 

P. Colangelo et al., Phys. Lett. B317 (1993) 183; S. Narison, Phys. Lett. 

B327 (1994) 354; J. M. Soares, Phys. Rev. D49 (1994) 283. 

A. Khodzhamirian, G. Stoll, and D. Wyler, Phys. Lett. B358 (1995) 
129. 

A. Ah and V.M. Braun, Phys. Lett. B359 (1995) 223. 

I.I. Balitsky, V.M. Braun, and A.V. Kolesnichenko, Nucl. Phys. 312 
(1989) 509. 

T.E. Browder and K. Honscheid, Prog. Part. Nucl. Phys. 35 (1995) 81. 

J.F. Donoghue, E. Golowich and A. A. Petrov, preprint UMHEP-433 
(1996) ||hep-ph/ 96095301 . 



[28] A. Ah, G. F. Giudice and T. Manuel, Z. Phys. C67 (1995) 417. 



21 



[29] S. Bertolini, F. Borzumati and A. Masiero, Phys. Rev. Lett. 59 (1987) 
180; R. Grigjanis et al., Phys. Lett. B213 (1988) 355; 

B. Grinstein, R. Springer, and M.B. Wise, Phys. Lett. 202 (1988) 138; 
Nucl Phys. B339 (1990) 269; G. Cella et al., Phys. Lett. B248 (1990) 
181. 

[30] M. Misiak in ref. §; 

A.J. Buras and M. Miinz, Phys. Rev. D52 (1995) 186. 

A. F. Falk, M. Luke and M. J. Savage, Phys. Rev. D49 (1994) 3367. 

A. Ah, G. HiUer, L.T. Handoko, and T. Morozumi, preprint DESY 96- 
206; Hiroshima univ. report HUPD-9615 ||hep-ph/96091ig| . 

A. Ah and E. Pietarinen, NucL Phys. B154 (1979) 519; G. AltareUi et 
al, NucL Phys. B208 (1982) 365. 

K. Fujikawa and A. Yamada, Phys. Rev. D49 (1994) 5890; 
P. Cho and M. Misiak, Phys. Rev. D49 (1994) 5894. 

A. Ah, T. Manuel and T. Morozumi, Phys. Lett. B273 (1991) 505. 

J. Hewett, Phys. Rev. D53 (1996) 4964. 

F. Kriiger and L.M. Sehgal, Phys. Lett. B380 (1996) 199. 

P. Cho, M. Misiak, and D. Wyler, Phys. Rev. D54 (1996) 3329. 

C. Albajar et al. (UAl), Phys. Lett. B262 (1991) 163. 

T. Skwarnicki, in Proc. of the 17th Int. Symp. on Lepton Photon In- 
teractions, Beijing, P.R. China, 10 - 15 August 1995; Editors: Zheng 
Zhi-Peng and Chen He-Sheng. 

CDF Collaboration, Fermilab Conf. 95/201-E(1995). 

A. Ah, C. Greub and T. Manuel, DESY Report 93-016 (1993), and in 
B -Physics Working Group Report, ECFA Workshop on a European B- 
Meson Factory, ECFA 93/151, DESY 93-053 (1993); Editors: R. Alek- 
san and A. Ali. 



22 



[43] G. Buchalla, A.J. Buras, and M.E. Lautenbacher, MPI-Ph/95-104; 

TUM-T31-100/95; FERMILAB-PUB-95/305-T; SLAG-PUB 7009; 
I ph/95l23"gl| . 

[44] Y. Grossman, Z. Ligeti, and E. Nardi, Nucl. Phys. B465 (1996) 369. 

[45] G. Buchalla, and A.J. Buras, Nucl. Phys. B400 (1993) 225. 

[46] Gontributed paper by the ALEPH collaboration to the International 
Gonference on High Energy Physics, Warsaw, IGHEP96 PAlO-019 
(1996). 

[47] GDF Gollaboration, Fermilab Gonf. 95/201-E(1995); P. Sphicas (GDF), 
private communication. 

[48] G.-L. Lin, J. Liu, and Y.-P. Yao, Phys. Rev. D42 (1990) 2314; H. 
Simma and D. Wyler, Nucl. Phys. B344 (1990) 283; S. Herrlich and 
J. Kalinowski, Nucl. Phys. 381 (1992) 501; T.M. Aliev and G. Tu- 
ran, Phys. Rev. D48 (1993) 1176; G.G. Devidze, G.R. Jibuti, and A.G. 
Liparteliani, preprint HEP Institute, Univ. of Tbilisi (revised version) 
(1995). 

[49] S. Abachi et al. (DO Gollaboration), contributed paper to the Interna- 
tional Gonference on High Energy Physics, Warsaw, IGHEP96 (1996). 

[50] L3 Gollaboration, contributed paper to the EPS Gonference, (EPS-0093- 
2) Brussels, 1995. 



23 



Decay Modes 


i3(SM) 


Measurements and 90% CL. Upper Limits 




(3.2 ±0.58) X 10-4 


(2.32 ±0.67) X 10-4 


{B^, B^) K*j 


(4.0 ±2.0) X 10-^ 


(4.2 ± 1.0) X 10-5 [| 


B') ^ X,7 


(1.6 ± 1.2) X 10-5 




5± ^ p± + 7 


(1.5 ± 1.1) X 10-*^ 


< 1.1 X 10-5 § 


5° ^ pO + 7 


(0.65 ±0.35) X 10-^ 


< 3.9 X 10-5 


i 


^ + 7 


(0.65 ±0.35) X 10-^ 


< 1.3 X 10-5 § 


fis ^ + 7 


(4.2 ±2.0) X 10-5 




Bs^K* + ^ 


(0.8 ±0.5) X 10-6 




(Bd, Bu) ^ X,e+e- 


(8.4 ± 2.2) X 10-6 




(5,,5„)^X,e+e- 


(4.9 ±2.9) X 10-^ 




^ X,p+p- 


(5.7 ± 1.2) X 10-6 


< 3.6 X 10-5 \m] 


{Bd, Bu) Xdfi+n' 


(3.3 ± 1.9) X 10"^ 




{Bd,B^)^Xsr+r- 


(2.6 ±0.5) X 10-^ 




{Bd,B^)^XdT+T- 


(1.5 ±0.8) X 10-^ 




iBd,B^)^Ke+e- 


(5.9 ±2.3) X 10-^ 


< 1.2 X 10-5 


m 


(Bd, Bu) ^ 


(4.0 ± 1.5) X 10-^ 


< 0.9 X 10-5 


m 


[Bd, Bu) ^ 


(4.0 ± 1.5) X 10-'^ 


< 0.9 X 10-5 


m 


{Bd,Bu)^ K*e+e- 


(2.3 ±0.9) X 10-6 


< 1.6 X 10-5 


m 


{Bd, ^ 


(1.5 ±0.6) X 10-6 


< 2.5 X 10-5 


m 


(Bd, Bu) Xs vv 


(4.0 ± 1.0) X 10-5 


< 7.7 X 10-4 


g 


{Bd, Bu) Xd vv 


(2.3 ± 1.5) X 10-6 




{Bd,B^)^Kvv 


(3.2 ± 1.6) X 10-6 




{Bd, Bu) ^ K* vv 


(1.1 ±0.55) X 10-5 




Bs 77 


(3.0 ± 1.0) X lO-'*" 


< 1.1 X 10-4 |50] 


Bd 77 


(1.2 ±0.8) X 10-^ 


< 3.8 X 10-5 


Bs ^ r+r- 


(7.4 ± 1.9) X 10-^ 




Bd ^ r+r- 


(3.1 ± 1.9) X 10-^ 




fi^ ^ P+;U" 


(3.5 ± 1.0) X 10-9 


< 8.4 X 10-6 


fid ^ p+p" 


(1.5 ±0.9) X 10-1° 


< 1.6 X 10-6 ig 


5, ^ e+e" 


(8.0 ±3.5) X 10-14 




fid ^ e+e- 


(3.4 ±2.3) X 10-15 





Table 1: Estimates of the branching fractions for FCNC -B-decays in the 
standard model taking into account the uncertainties in the input parameters 
as discussed in . The entries in the second column correspond to the short- 
distance contributions only, except fo?4he radiative decays B^" p^±7 and 
B^ (p°, C(j)±7, where long-distance effects have also been included. For the 
two-body branching ratios, we have used fs^ = 200 MeV and fss/ fsa = 1-16. 
Experimental measurements and upper limits are also listed. In the second 
row, the statistical and systematic uncertainties have been combined to give 
the quoted experimental uncertainty. 



